Realizing non-Abelian gauge potentials in optical square lattices: 
Application to atomic Haldane-Chern insulators 
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We describe a scheme to engineer non-Abelian gauge potentials on a square optical lattice using laser-induced 
transitions. We emphasize the case of two-electron atoms, where the electronic ground state is laser coupled 
to a metastable state in a state-dependent optical lattice, and discuss how the nuclear spin can be viewed as a 
"flavor" quantum number undergoing non-Abelian tunneling along nearest-neighbor links. We show that this 
technique can be useful to simulate the equivalent of the Haldane quantum Hall model using cold atoms trapped 
in a square optical lattice, offering an interesting route to realize Chern insulators. The emblematic Haldane 
model is particularly suited to investigate the physics of topological insulators, but requires, in its original form, 
complex hopping terms beyond nearest-neighboring sites. This drawback inhibits a direct realization with cold 
atoms, using standard laser-induced tunneling techniques. We demonstrate that a simple mapping allows to 
express this model in terms of matrix hopping operators, that are defined on a standard square lattice. This 
mapping is investigated for two models that lead to anomalous quantum Hall phases. We discuss the practical 
implementation of such models using laser-induced tunneling methods. 



I. INTRODUCTION 

The discovery of the quantum Hall (QH) effect, in two- 
dimensional electronic systems subjected to large magnetic 
fields, revealed the existence of novel quantum phases: the 
topological insulating states [1, 2]. Such phases are insulating 
in the bulk but exhibit current-carrying modes on the edge 
of the sample. These chiral propagating modes, whose 
energies are located within the bulk gaps, are protected by 
topology. They can be formally related to the existence 
of non-trivial topological invariants - Chern numbers - 
associated with the band structure [3, 4]. The existence of 
edge modes is guaranteed as long as the topology of the 
band structure does not change, irrespective of the details of 
the microscopic Hamiltonian. In the early age of quantum 
Hall physics, it was generally thought that this phenomenon 
could only emerge in two-dimensional systems featuring 
Landau levels, namely, in samples necessarily subjected to 
magnetic fields. Haldane showed in 1988 that this was not a 
necessary condition [5]. To this aim, he introduced a simple 
lattice model exhibiting the quantum Hall effect without 
Landau levels, thereafter referred to as the anomalous QH 
effect. In his seminal work [5], Haldane showed that quantum 
Hall phases were rooted in the breaking of time-reversal 
symmetry, which potentially increased the possibilities to 
access this phenomenon in a wider range of two-dimensional 
systems. This first step played a major role in the quest 
for new topological states, as it paved the way for the 
generalization of the quantum Hall effect. Indeed, this work 
inspired studies aiming to discover topologically-ordered 
band structures in other spatial dimensions, and also, in 
physical systems satisfying time-reversal symmetry [6]. This 
effort led to the prediction [7-9] and to the discovery [10, 11] 
of Z 2 topological insulators, which are materials featuring 
large spin-orbit coupling and exhibiting the quantum spin 
Hall effect. In fact, the transport properties offered by these 



materials, which currently attract the scientific community for 
their potential future applications, can be described by simple 
lattice models that are direct generalizations of the initial 
Haldane model [7, 12]. Besides, it was recently shown that 
the Haldane model could be tailored in order to create band 
structures presenting low-energy flat bands with non-zero 
Chern numbers [13]. In such configurations, in analogy with 
fractional quantum Hall states occurring in strong magnetic 
fields, interactions can profoundly change the many-body 
ground state. This offers a novel route towards strongly- 
correlated topological liquids, the so-called fractional Chern 
insulators [13-16]. Although it played a fundamental role 
in the search for topological materials, the Haldane model 
was never reproduced in solid-state systems. Recently, it was 
suggested that the anomalous QH effect could be realized 
and observed in HgMnTe quantum wells [17], in silicene 
subjected to in-plane magnetic fields [18], and in cold-atom 
setups [19-26]. 



The original Haldane model [5] is defined on a honeycomb 
lattice with real nearest-neighbor (NN) and complex next- 
nearest-neighbor (NNN) hoppings, see Fig. 1 (a). The NN 
hopping alone leads to a gapless energy spectrum featuring 
two Dirac cones familiar from graphene studies [27]. The 
complex NNN hopping, represented by red and blue arrows 
in Fig. 1 (a), opens a single gap at half-filling (i.e. around 
E = 0). This bulk energy gap is non-trivial [5]: when the 
Fermi energy Ep w lies in this gap, a single edge mode 
is populated and carries current around the system [4], lead- 
ing to a quantized Hall conductivity, anih/e 2 ) = v = ±1. 
From a topological point of view, the lowest energy bulk band 
E(k) < is associated with a non-vanishing Chern number 
v = ±1, which guarantees the presence of the robust and 
unique edge mode within the bulk gap [3, 4, 28]. The Chern 
number v can be directly evaluated from the effective Dirac 
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equations associated with the two Dirac cones [5, 27], 

v = i(sign(M K +) - sign(M K -)), (1) 

where M K ± are the effective masses related to the two 
independent Dirac points K*. 

Realizing the Haldane model with cold atoms trapped in 
an optical lattice is attractive, as it would provide a simple 
playground to investigate the physics of topological (Chern) 
insulators. Such a proposal was described in Refs. [23, 24], 
where laser-induced tunneling [29] was considered to couple 
two state-dependent triangular optical lattices. The simplicity 
of the Haldane model relies on the fact that it is a two-band 
model, featuring a unique energy bulk gap and hosting a 
single topological edge mode (y = ±1). This minimal 
topologically-ordered system is to be compared, for instance, 
to the emblematic Hofstadter model [30], i.e. a tight-binding 
model for an electron in a uniform magnetic field moving on a 
square lattice, which is characterized by a multi-band energy 
spectrum and hosts (infinitely many) quantum Hall phases 
with arbitrary v <G Z [31]. In this sense, the two-band Haldane 
model is mathematically easier to handle than the Hofstadter 
model [24], and it thus constitutes a good basis for studying 
the effects of interactions, in particular, in view of realizing 
fractional Chern insulators [13-16, 32] and topological Mott 
insulators [26] with cold atoms. Finally, the Haldane model 
constitutes the building blocks for the Kane-Mele model [7], 
and therefore, its optical-lattice realization opens a possible 
route for the observation of the quantum spin Hall effect with 
cold atoms [33-36]. 

In this work, one exploits the fact that the honeycomb 
Haldane model can be mapped into a two-component square 
lattice with NN and NNN matrix hoppings. We demonstrate 
that the NNN hoppings, which are inconvenient for an optical 
lattice implementation, can be simply omitted: they do not 
contribute to the appearance of a non-trivial topological 
phase. This simple observation makes the realization of 
the Haldane model feasible with laser-assisted tunneling 
on a square lattice, at the price of an increased complexity: 
the resulting model involves non-Abelian gauge potentials 
[37-39]. The models stemming from this non-Abelian 
framework are versatile, and thus, they can be generalized to 
simulate an assortment of topological systems and Dirac-like 
systems [33,40,41]. 

In Section II, we first describe an extension of the scheme 
proposed in [37] allowing one to realize such non-Abelian 
gauge potentials. In Section III, we show how these synthetic 
gauge potentials can be tailored in order to reproduce the orig- 
inal Haldane model on a square lattice. We also briefly com- 
ment on the detection of relevant signatures, based on avail- 
able experimental probes. In Section IV, we show that similar 
considerations can be applied to the so-called 7r— flux model 
[13, 38, 42], which leads to identical physics with a simpler 
setup. 



II. NON-ABELIAN GAUGE POTENTIALS FOR 
ULTRACOLD ATOMS IN OPTICAL LATTICES 

We start by outlining an experimental scheme generating 
effective non-Abelian gauge fields for ultracold fermionic 
atoms on a lattice. The scheme is a generalization of the ones 
proposed by Osterloh et al. [37], building on the earlier pro- 
posal by Jaksch and Zoller to realize effective magnetic fields 
[29]. Both proposals [29, 37] are designed for alkali atoms, 
and rely on using multilevel atoms in a spin-dependent opti- 
cal lattice, i.e. a periodic trapping potential made from several 
different sublattices trapping different internal states. A set 
of laser beams drive Raman transitions that change the inter- 
nal state, thereby inducing tunneling from a given sublattice 
to another. Crucially, the laser-assisted tunneling process is 
characterized by a complex matrix element, with an argument 
determined by the laser phase. A suitable configuration of 
laser beams then leads to a non-vanishing geometrical phase 
for an atom tunneling around a closed trajectory, which real- 
izes an effective magnetic field. 

Implementing this scheme experimentally with fermionic 
alkali atoms is unfortunately plagued by spontaneous emis- 
sion: the spin-dependent lattice and Raman lasers frequencies 
need to be set relatively close to resonance, leading to heating 
and losses hardly compatible with quantum gases experiments 
[43]. Moreover, such spin-dependent lattices are sensitive to 
the magnetic moment of the atoms, leaving the protocol vul- 
nerable to stray magnetic fields. These issues can be avoided 
by using alkaline earth atoms (or atoms with a similar level 
structure, such as Ytterbium). We outline below how the pro- 
posals must be modified to use these atoms. We first start 
by recalling the scheme proposed in Ref. [43] to realize an 
effective magnetic field - or equivalently an Abelian gauge 
structure with U(l) symmetry - and we explicitly apply it to 
a checkerboard geometry, see Fig. 2 (a). We then explain in 
a second step how this can be extended to non-Abelian con- 
figurations in the spirit of Osterloh et al. [37]. For simplicity 
we consider an atom with spin 1/2 and generations of SU(2) 
gauge potentials. In principle, the method can be extended 
to other gauge groups but this typically requires many more 
lasers than the simpler U(l) or SU(2) schemes we outline here. 
As a result, the experimental implementation becomes signif- 
icantly more involved. 



A. Spin-dependent optical lattices, laser-induced tunneling 
and the Peierls phase 

Two-electron atoms (such as alkaline earth atoms or Ytter- 
bium) feature two peripheral electrons that form a singlet 1 Sq 
ground state and several low-lying triplet states, including a 
long-lived metastable 3 P state (with a typical lifetime of tens 
of seconds). From now on, the notation g =1 5 and e = 3 Pq 
will be used for simplicity. The ultra-narrow g — e transition 
enables a direct coherent manipulation of the internal states 
without spontaneous emission. 

Consider a two-dimensional square lattice populated with 
atoms in the two internal states g and e. Atoms in both states 
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(a) Haldane model (b) Haldane model with matrix hopping 




Figure 1: (a) The Haldane model on the honeycomb lattice. The unit cells of the honeycomb are labeled by the coordinates (m,n) G Z. 
Two inequivalent sites A, B belonging to the same unit cell are connected by a full black line. Standard nearest-neighbour (NN) hopping, 
with amplitude J and denoted by dotted lines, takes place between the A (red) and B (blue) sites. The complex next-nearest-neighbor (NNN) 
hoppings, introduced to open a quantum Hall gap, are represented by thick blue and red lines. The NNN tunneling factors are +i\ according 
to the orientation designated by the circular arrow, and — iX otherwise, i.e., they introduce a chirality in the system, (b) The same Haldane 
model translated into a non-Abelian square lattice, with matrix hopping operators U x ^ y i. The "undesired" diagonal hoppings D are depicted 
by red dotted arrows and disappear in the limit a = (cf. text). Note that when a — 0, this model reduces to the non-Abelian optical lattice 
illustrated in Fig. 2 (c). 



(a) State-dependent lattice (b) States and transitions (c) Non-Abelian optical lattice 




Figure 2: (a) Sketch of the spin-dependent trapping potential for atoms with nuclear spin 1/2. Atom in electronic state g, shown as green 
circles, (resp. e, shown as dark red circles) are trapped at the potential minima (resp. maxima) independently of their nuclear spin. This results 
in a spatial arrangement on a checkerboard pattern. A laser resonant on the g — e transition can induce tunneling along the diagonals e^, 
By. We show for illustration a closed trajectory around a unit cell of the checkerboard lattice, (b) An applied magnetic field shifts the various 
transitions between internal states depending on the value of the nuclear spin, allowing independent addressing of each of them (for instance, 
addressing the rightmost it transition independently from the a~ transition), (c) The resulting "non-Abelian" optical lattice, with U(2) hopping 
operators U x i_ y i qt 12x2 acting along nearest-neighboring sites. Note that we use the notation (m,n) [resp. (m, n)] to designate the lattice 
sites in the x' — y [resp. x — y] axis system. 
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are confined to the x — y plane, independently of the inter- 
nal state, by a laser set near a "magic" wavelength where both 
states have equal polarizabilities. Along the x and y direc- 
tions, two independent standing waves are applied to form a 
square two-dimensional optical lattice. A state-dependent op- 
tical lattice is realized if the trapping laser is set at an "anti- 
magic" wavelength, such that the polarizability for states g 
and e are opposite: negative for g atoms which are trapped 
near the intensity maxima, positive for e atoms which are 
trapped near the minima. This results in a lattice poten- 
tial with checkerboard geometry, where sites corresponding 
to the e sublattice are centered on the cells of the g sublat- 
tice (see Fig.2 (a)). In the following it is convenient to con- 
sider a rotated frame, with axes e' y = (e x + e y ) /\[2 and 
e' x = (e x — e y ) /\pl corresponding to the diagonals of the 
original square sublattices. 

Tunneling between g/e sublattices, with amplitude J, is 
enforced by a near-resonant laser beam propagating along y, 
with Rabi frequency ft and wavevector q = qe y . As shown 
below, the coupling matrix element hQe lqv realizes the so- 
called "Peierls-Luttinger substitution" J — > Je lA , describ- 
ing how a gauge potential A modifies the particle tunneling 
J (and hence, the band structure) in the tight-binding limit 
[30, 44]. Following Jaksch and Zoller [29], we assume the 
atoms only occupy the fundamental Bloch band and make 
the tight-binding approximation. We use a shorthand nota- 
tion (m, n) for the site located at x = md/2, y = nd/2, where 
m, n are even for state g and odd for state e, and where d is 
the lattice spacing for each sublattice, see Fig. 2(a). We then 
obtain the effective g — e hopping matrix element between two 
neighboring sites, located at r g = (m, n) and r e = r g + e' x , 
as 

J(x) e iq-r. = ^ e i™ J w * {r _ e > x)e i*r Wg{r) ^ (2) 

with Ct and q = qe y the laser Rabi frequency and wavevector, 
respectively. The phase factor has been written as 2na = qd, 
where a can be controlled from zero to a value larger than 
one by modifying the orientation of the coupling laser. The 
magnitude of the tunneling matrix element is controlled by 
the Rabi frequency and the overlap integral between the 
Wannier functions w g / e associated with each sublattice. 
In principle, the laser is also able to induce tunneling to 
next-nearest-neighbors of the composite lattice, and beyond. 
In practice, the matrix elements are exponentially suppressed 
compared to the ones describing nearest-neighbors tunneling, 
and negligible for realistic experimental configurations. This 
sets a constraint on the class of models that can be realized in 
this way, as we will discuss later in the paper. 

Let us note that the checkerboard lattice considered here po- 
tentially allows to induce the tunneling and the Peierls phases 
along both spatial directions (in contrast with the "column" 
geometry previously considered in Refs. [29, 43, 45]). This 
central ingredient will be largely exploited in the following of 
the present proposal, see Section II C. 



B. Realizing Abelian gauge potentials 

First of all, we note that a particle making a loop around a 
unit cell of the checkerboard lattice with tunneling elements 
(2) does not acquire any phase factor. This is a consequence 
of unitarity, which imposes that the laser coupling matrix el- 
ements for g — > e transitions [see Eq. (2)] and its reversed 
e — > g counterpart, are complex conjugates, namely 

r g ^r e = r g +e^. J#> e*-- , 

r e = r g +e'^r g : (j£>)* e~^' . (» = x, y) 

(3) 

Because of the alternating pattern of g —¥ e and e — > g links 
[see Fig. 2 (a)], the sign of the tunneling phases alternates 
when moving in a given direction, which indeed leads to the 
cancellation of the overall phase picked by a particle hopping 
around a unit cell of the checkerboard. Starting from the site 
(m, n), and using Eqs. (2)-(3), one obtains that the phase ac- 
quired along the closed path (m, n) — > (m + 1, n + 1) — > 
(m + 2, n) -> (m + 1, n - 1) -> (m, n) is 

27t$q = nan — nan + nan — nan = 0, (4) 

which indicates that the effective "magnetic" flux <£>□ pene- 
trating each plaquette (in units of the flux quantum) is zero. 
This trivial flux configuration is in contrast with the "column" 
geometry discussed in Refs. [29, 43, 45], where a staggered 
effective magnetic flux is generated, and it results from the 
higher symmetry of the checkerboard lattice. This can be 
cured by an additional superlattice with double period 2d, gen- 
erating a potential of the form V"(x) oc cos 2 (nx/4:d + ip). A 
suitable choice of the relative phase ip leads to on-site energies 

E g - AV, for x/d = 4m, 
E e , for x/d = 4m + 1, 
Eg, for x/d = 4m + 2, 
E e + AV, for x/d = 4m + 3, 

with m an integer. Here the E a (a = e,g) denote the in- 
ternal energies in free space, and AV is the amplitude of the 
modulation with period 2d x . In this potential landscape, the 
resonance frequencies for transitions linking neighboring sites 
become non-degenerate, 

wi = uj q + AVyftfor (4m, 4n) (4m + l,4n + 1), 
lo 2 = W!- AV/h for (4m + 1, 4n + 1) ->• (4m + 2, 4n), 
u> 3 = u>! for(4m + 2,4n) -> (4m + 3,4n + l), 

uj 4 = ui + AVyftfor (4m + 3,4n + 1) ->• (4m + 4,4n), 

with hbJo = E e — Eg the bare transition frequency. 

To the state-dependent lattice, one thus applies three cou- 
pling lasers propagating along y. The laser at frequency uii 
is chosen with a wavevector q = qe y , and the lasers at fre- 
quencies u)\ ± AV with the opposite wavevector q' = —qe y . 
If we neglect off-resonant transitions, the alternation of the 
wavevectors compensates the alternation of the sign of the 
tunneling phases, thereby leading to a non-zero "magnetic" 
flux per plaquette that is uniform across the lattice. 
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C. Realizing non-Abelian gauge potentials 

This proposal can be generalized to non-Abelian gauge po- 
tentials when 5 and e now represent manifolds of degenerate 
states (see Fig. 2 (a)-(b)). Let us consider a 2 x 2 system with 
states {5-1/2,5+1/2} and {e_ 1/2 , e +1/2 }, as realized in the 
fermionic 171 Yb with nuclear spin 1/2. In a moderate mag- 
netic field (a few tens of Gauss), the Zeeman effect allows 
to distinguish between <r T transitions 5+1/2 — > e =pi/2, that 
change the spin projection, from tt transitions 5+1/2 — > e±i/ 2 
that conserve it (see Fig. 2 (b)). This allows to correlate tun- 
neling in a spatial direction with rotations in internal states 
and state-dependent tunneling phases. Laser-assisted tunnel- 
ing along x' or y' is now described by 2 x 2 matrices U x , 
and U y i acting on a two-component spinor. If it is possible to 
arrange such that [U x ,, U y >] ^ 0, this mimics a non-Abelian 
gauge potential A through U = e*f A ' dl [57]. 

In order to realize U x , ^ U y >, one must distinguish between 
the tunnelings directed along e' x and e' y . Therefore, this re- 
quires two superlattice potentials along the x and y axis, with 
different heights AV X and AV y . This generates nine differ- 
ent transition frequencies differing by 0, AV x ,AV y , ±AV X ± 
AV x ,2AV x ,2AV y , that can be grouped into two families. 
One family corresponds to transitions along e' x and the other 
one corresponds to transitions along e' y . 

We first consider the transitions along e' y for illustration. 
Let us choose the coupling lasers propagating along z (thus 
inducing tunneling with no relative phase for motion in the x— 
y plane), with a linear polarization a x , and tuned on resonance 
in the absence of magnetic field. Such lasers will induce a 
tunneling matrix oc a x . A second set of lasers is added, with 
polarization <r v , on resonance with transitions along e' x . This 
configuration generates a spin-orbit coupling of the form 

U x , oc & y , U y , oc e^&x, (5) 

where the relative phase <j>i of the two sets of lasers determine 
the relative sign of the tunneling matrices. 

More complicated tunneling matrices with complex 
position-dependent elements can be generated using a differ- 
ent laser configuration. Two specific configurations, leading 
to the analog of the Haldane model [5], will be discussed be- 
low. In principle, arbitrary tunneling matrices can be gener- 
ated in this way. However, the increase in complexity makes it 
experimentally very challenging. The configuration described 
above, for instance, generates no less than 9 different frequen- 
cies. Those frequencies are typically in the 10 kHz range, 
much less than the absolute optical frequency uj correspond- 
ing to the bare atomic transition. As a result, they can be 
generated from a single laser using acoustic-optical or electro- 
optical modulators. This represents nevertheless a significant 
technological challenge for designing arbitrary tunneling ma- 
trices. 



III. THE HALDANE MODEL ON THE SQUARE LATTICE 

A. Mapping the Haldane model to a non-Abelian model on a 
square lattice 

The idea that we expose in this Section is to propose and 
analyze a mapping from the original Haldane model, involv- 
ing complex tunneling phases on a honeycomb lattice [5], into 
a spin-1/2 model defined on the square lattice, see Figs. 1 (a)- 
(b). This is possible in two dimensions, since any lattice with 
N inequivalent sites within its unit cell can be formally la- 
belled using: (a) two spatial coordinates (to, n) describing the 
location of the unit cell, and (b) an additional "pseudo-spin" 
index r = 1, . . . , N labeling the inequivalent sites. Since the 
honeycomb lattice features two inequivalent sites, denoted A 
and B, it can be mapped into a square lattice hosting spin-1/2 
objects. 

Our starting point is to describe the honeycomb lattice 
using the notations illustrated in Fig. 1 (a). In these notations 
[42], one uses discrete coordinates (to, n) to label the unit 
cells of the honeycomb, where to, n are integers. Here 
each "site" (to, n) hosts the wave-functions ^4 (to, n) and 
i/j B (m, n), associated with the two inequivalent sites of the 
honeycomb. In other words, the honeycomb lattice can 
be interpreted as a square lattice with a pseudo-spin 1/2 
structure. Thus, the tight-binding Hamiltonian describing 
NN hopping between A and B neighboring sites can be 
mapped into an equivalent operator, which involves hopping 
between opposite pseudo-spin components, which are located 
at neighboring sites of a square lattice, see Fig. 1 (b). This 
mapping smoothly deforms the original Brillouin zone and 
energy bands but does not affect the physical properties, such 
as bulk gaps and topological phases, which are discussed in 
the following. 

One considers the Haldane model, which features both NN 
and NNN hoppings on the honeycomb lattice (Fig. 1 (a)). The 
single-particle Schrodinger equations, satisfied by the wave- 
functions Va(w, n) and ip B {m, n), take a very simple form 
in terms of the notations introduced above, and read 

EipA(m,n) = ip B {m,n) + ij) B {m- l,n) + ipB(m, n - 1) 

+ iX (^4>a {m + 1, n) — aipA{m + 1, n — 1) + ipA^m, n—1) 

- ^a{tti - 1, n) + ai/j A (m - l,n + 1) - i/) A (m,n + 1)^ , 
Eip B (m,n) = ip A (m,n) +^ A {m + l,n) +ip A (m, n+ 1) 
+ iX ^V'B( m , n + 1) — V'-B( m + 1) n ) + aipBim + 1, n — 1) 

- tp B {m,n - 1) + tpB{m - l,n) - aip B (m - l,n + l)j , 

(6) 

where the tunneling amplitude J to nearest-neighbors of the 
honeycomb lattice is set to J = —1, and thus defines our unit 
of energy. Here, A is the amplitude of the NNN hoppings on 
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the honeycomb lattice and it corresponds to the ti parameter 
of the original Haldane model (or equivalently to the intrinsic 
spin-orbit coupling strength in the spin- 1/2 Kane-Mele model 
[7]). We remind that the complex NNN hoppings, which 
are responsible for the opening of the non-trivial topological 
bulk gap in this model, introduce a chirality in the system, 
see Fig. 1 (a). Note that we have included an additional 
parameter a e [0,1]- for reasons that will appear clear after 
the mapping to a square lattice - with a = 1 in the original 
Haldane model. 

Setting ^ (m, n) = (^(m, n), V>b(?ti, n)), one can rewrite 
Eq. (6) as 

£7* (to, n) = F^(m, n) + U x ^(m + 1, n) + U y ^(m, n + 1) 
+ £>*(m + l,n- 1) + Ul^(m - l,n) 
+ Ul,*(m,n-l) + D^(m-l,n+l), (7) 

where U x > = b- + i\<r z , U y > = <r_ — i\a z , F = a x and 
D = —ia\a z = —it$ m% b z . Here, we have introduced the 
Pauli matrices a x ,y,z an d = i?x ^ i&y)/2- Equation (7) 
describes a non-interacting two-component system, evolving 
on a square lattice, see Fig. 1 (b). The obtained model is 
characterized by non-Abelian hopping operators U x ^ y ' that 
act along the links, but also by an onsite spin-mixing term 
F and diagonal hoppings D with amplitude t^. ig = c*A, see 
Fig. 1 (b). To keep with the notations used in the preceding 
Section, we denote by x' , y' the principal axes of the square 
lattice. 

Realizing non-Abelian hoppings U x i tV i between the 
nearest-neighbors of a square optical lattice is a possible but 
difficult task, as discussed above. The most important issue 
is the presence of additional diagonal matrix hoppings D, 
which will be exponentially small in a cold atom realization 
compared to f/ x /y (see discussion in the preceding Section). 
Consequently, a natural question arises: how important is 
this diagonal hopping D for the obtention of the topological 
phase? Would the non-trivial bulk gap and corresponding 
edge-states survive in the limit a — > 0? We answer this 
question by investigating the fate of the energy spectrum 
and topological order (Chern number and edge-states) as the 
Haldane model described by Eq. (7) is reduced to the simpler 
model with D = 0. 

In momentum space, the Hamiltonian describing our sys- 
tem reads 

H{k x ,,k y ,) = d{k)-a, (8) 

where a is a vector of Pauli matrices and where the vector 
d(fc) is given by 

(cos(£v) + cos(Ay) + 1 \ 
sin(Av) + sin(Ay ) 
2A [sin(Ay) — sin(Av)] + 2aXsm(k x i — k y i) J 



In the limit A = 0, the energy spectrum is gapless and de- 
scribes a semi-metal at half-filling. The two branches are 

E ± (k; \ = 0) = ± | d(fc; A = 0) |, (9) 

= ±y^3 + 2(cosfc X ' + cos k y i + cos(k x i — Ay)), 

and they touch at the two inequivalent Dirac points K + = 

(27r/3,47r/3), K = (47r/3,27r/3), where E± = (or 
equivalently, where d x + id y = 0). When A/0, a gap 

A = 2|4(if + )|-4|A|^(1 + |) (10) 

opens at the Dirac points. It is well established for two- 
component lattice models that the topology of the ground band 
is entirely characterized by the vector d parameterizing the 
Hamiltonian: the Chern number of the band v is identical to 
the winding number of d(fc), as k is varied in the first Bril- 
louin zone [12, 23, 24, 28]. As already presented in Eq. (1), 
this winding number v can be directly related to the effective 
masses associated with the two Dirac points K ± . Here, these 
masses are given by 

M ± = d z (K ± ) = T 2XV3(l + ^), (11) 

as can be deduced by developing the Hamiltonian H(k x > , Ay ) 
in the vicinity of the Dirac points K ± . The effective masses 
have different signs, leading to a non-zero Chern number 
v = ±1. Moreover these masses are non-vanishing in the 
limit a — > 0. Therefore the bulk energy gap opened by the 
perturbation A, stemming from the complex Haldane NNN 
hopping, survives when a = 0. Since the masses pre- 
serve their sign during this transformation, the Chern number 
in Eq. (1) remains constant and non-trivial. The latter remark 
can also be formulated in the following manner: since the bulk 
gap is preserved during the transformation a = 1 — > a = 0, 
the topological Chern number characterizing the bulk bands 
is unaffected. This result is further illustrated in Fig. 3, 
which shows the energy spectrum for the two limiting cases 
a = and a = 1. In Figs. 3 (a)-(b), the spectrum was 
obtained using a cylindrical geometry aligned along x'\ it 
shows the projected bulk bands E±(k) — > E±(k y ') and 
the topological edge-states inside the bulk gap [4]. This 
result demonstrates that the diagonal hopping induced by the 
operator D = —ia\a z in Eq. (7) does not play any role in 
the realization of topological insulating phases and can thus 
be omitted for the sake of experimental feasibility. 

To summarize, the physics of the Haldane model can be 
explored with an optical square lattice described by the sim- 
plified tight-binding Hamiltonian 

^ = ^y^j C ln+l,rftx' c m,n + c ln,n+lUy' c m,n + c L-l,n^x' C ™*>™ 
C m,n— lUy/Cm.n + c m.n^ c in,n! (12) 

U x i = <7_ + iAcTz, U y i = cj_ — i\(7 z , F = a x , 

where the two-component operator c\ n n creates a particle at 
site (m, n) and where the spin- 1/2 structure is implicit. We 
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Figure 3: Energy spectrum E = E(k y i) for A = 0.05t: (a) a = 1, 
(b) a = for a cylindrical geometry aligned along x'. This spectrum 
shows the projected bulk bands E± (k) — > E± (k y i), and also reveals 
the presence of topological edge-states inside the bulk gap [4], In the 
standard case a = 1, the topological bulk gap is A = 6\/3A « 0.5. 
In the absence of the diagonal matrix hopping, i.e. a = 0, the gap 
is A = 4\/3A « 0.34 but survives. The energy unit is given by the 
tunneling amplitude J. 



where we used Eq. (2), and where n designates the y 
coordinate of the site (to, n), which is here labelled according 
to the x' — y' axis system, see Fig. 2. We note that each phase 
4>_ (to, n) is driven by a specific laser, corresponding to the 
different transition frequencies offered by the superlattice 
potential. By adjusting the relative phases of these lasers, 
one can annihilate the undesired phases <j>_ (to, n) = 
uniformly on the lattice, and thus realize constant tunneling 
operators U X ' ,y' oc <r_ as in Eq. (12). Combining the lasers 
generating the tt and er~ transitions, one obtains the required 
NN tunneling terms in Eq. (12). Finally, the on-site term 
F oc a x can be generated by an additional radio frequency 
field resonant at the Larmor frequency, acting uniformly on 
all sites [58]. 

Although the scheme seems directly feasible, it is rather 
complex due to the large number of transition frequencies in- 
volved. In the following Section IV, we present an equivalent 
model that leads to a simpler implementation. 



note that this system is formally similar to the HgMnTe quan- 
tum wells model proposed in [12, 17] to realize the anomalous 
quantum Hall effect. 

B. Implementation 

We now discuss how the Hamiltonian Eq. (12) can be im- 
plemented using the methods discussed in Section II. We start 
by the terms U x > — +i\a z and U' , = —i\a z . First, we con- 
sider transitions corresponding to tunneling along the x' direc- 
tion. For each value of the nuclear spin (+1/2 or —1/2), one 
can apply laser-assisted tunneling with tt polarization. The 
lasers propagate along z (thus inducing no tunneling phase), 
with a polarization along y parallel to the direction of an ap- 
plied magnetic field large enough to split the different transi- 
tions. Since the transition frequencies for g 1 / 2 — > and 
<7_i/2 — > e_i/ 2 are different, they must be addressed sepa- 
rately with independent lasers, see Fig. 2 (b). We choose a 
phase of ±7r/2 for the two lasers driving the two tt transitions 
for nuclear spins ±1/2. This gives a tunneling matrix propor- 
tional to 

U x > oc °.) =ia x . 

Another set of lasers can be used to generate tunneling op- 
erators along e' y , also driving tt transitions for nuclear spins 
±7r/2 with phases =F7r/2. This generates a tunneling matrix 

Uyi oc — ia z . (13) 

Next, we consider additional coupling lasers propagating 
along +e y , with a~ polarization, thus generating a tunneling 
matrix oc <r_ . For instance, the tunneling matrix element for 
the transition linking site r g = (m, n) to r e = r g + e' x = 
(to + 1, n) is then 

(m + l,n\ U x , | to, n) oc e i<t>-(™,n) a _ = e i™<* & _^ (14) 



C. Detection 

The Haldane-like optical lattice described above is charac- 
terized by the following physical properties: 

• a massless Dirac-like spectrum for A = 0, 

• a massive Dirac-like spectrum for small A< J, 

• a largely gapped spectrum for A <~ J, 

• Chern insulating phases for A ^ 0. 

The Dirac-like spectrum can be experimentally detected 
along the lines of the recent experiments [46, 47]. In par- 
ticular, the techniques developed in these experiments could 
identify the opening of the bulk gap in the vicinity of the 
Dirac points, and thus, evaluate the amplitude of the mass 
terms M^ 1 ^ 0. We point out that the detection of a Dirac 
spectrum in our square optical lattice setup would already 
be an interesting signature of the non-trivial synthetic gauge 
potential [38]. 

The large gap opened by the terms proportional to A could 
be directly detected through in-situ imaging. Indeed, this 
spectral gap directly manifests itself in the spatial density of 
the atomic cloud [43]. We show the corresponding density 
for A = and A = 1 in Fig. 4, for a system trapped by 
an additional harmonic potential V(r) = Vo(r/ro) 2 : the 
clear plateau depicted by the density when A = 1 is a direct 
signature of the gap opening. 

The most efficient way to identify the formation of a Chern 
insulating phase would be to detect the presence of robust and 
chiral edge modes, whose energies are located within the bulk 
gap discussed above. This could be performed by deforming 
the atomic cloud and imaging its time-evolving spatial density 
[48], or using light scattering methods [49, 50]. In principle, 
measuring the Chern number v could be performed through 



8 



(a) (b) 




Figure 4: (a) The spatial density for the Haldane-like optical lattice 
described by Eq. (12), for A = and A = 1. The atoms are trapped 
by a harmonic potential V(r) = Vo(r/ro) 2 , with Vb = 8J and 
ro — 17a, and the lattice size is 40 x 40. The density was obtained 
by setting the Fermi energy to the value Ep = 3 J (corresponding to 
the maximum of the highest bulk band for A = 0). In this figure the 
n coordinate is chosen to be fixed at the center of the trap, (b) Same 
as (a) for A = 1, but represented in the 2D plane x' — y . Note the 
clear plateau for A = 1, which indicates the opening of a spectral 
gap [43]. 



spin-resolved momentum density measurements [23, 24] or 
Bloch oscillations [51]. 



IV. THE 7T-FLUX MODEL 

The physics of the honeycomb Haldane model can be al- 
ternatively studied by considering non-trivial hoppings on the 
7r-flux model [13, 38, 42]. The latter is illustrated in Fig. 5 (a) 
and is characterized by tunneling amplitudes ± J along alter- 
nate columns. In the absence of additional hoppings, this sys- 
tem exhibits a gapless spectrum with two Dirac cones, and it is 
therefore physically equivalent to the honeycomb lattice [42]. 
In order to open a non-trivial topological gap, with Chern 
number v = ±1, complex diagonal hoppings are required, 
see the purple arrows in Fig. 5 (a). As for the original hon- 
eycomb model discussed in the previous Section, these chiral 
and complex NNN hoppings highly diminish the feasibility 
of this Haldane-like model with cold atoms trapped in optical 
lattices. Note that the 7r— flux model without NNN hoppings 
is equivalent to the Hofstadter model with the magnetic flux 
set to $= 1/2 [43, 52]. 

In this Section, we follow the same strategy as before: 
we consider the two-component wave-function ^>(m,n) — 
[ipA(n r i,n),'tljB(fn,n)), which describes particles on alter- 
nate columns, see the blue (A) and red (B) sites in Fig. 5 
(a). In these notations, the single-particle Schrodinger equa- 
tion satisfied by fy(m, n) in this model then takes the form 

E^(m, n) = F^>(m, n) + U x ^(m + l,n) + U y >^(m, n + 1) 

+ Di^(m + 1, n + 1) + D 2 #(m + 1, n - 1) 

+ Ul,^(m - l,n) + U^(m,n- 1) 

+ t)\^(m - 1, n - 1) + D\^(m - 1, n + 1), 

(15) 



where U x i = <r_, U y i = —<j z — i\a x , F — a x and 
L>i 2 = ±iaA<7_. Therefore, the model is translated into a 
non-Abelian square lattice, with direct tunneling operators 
U x ', y ', an onsite term F, and two diagonal hopping matrices 
£»i j2 , see Fig. 5 (b). 

The 7T— flux model can be analyzed in the same way as in 
the preceding Section, with identical conclusions. In the limit 
A = 0, the energy spectrum is given by 

E± (A = 0) = ±y / 2 + 2(cos£v +2cos 2 k y ,), (16) 

with two independent Dirac points at K + = (ir,ir/2) and 
K = (tt,3tt/2). In the vicinity of the Dirac points, the 
effective masses are given by 

M ± = ±2A(1 + a). (17) 

As a result, the Chern number v = ±1 remains non-zero even 
without diagonal hoppings (a = 0), as found previously. 

The 7T— flux quantum Hall model can thus be explored by 
realizing the tight-binding Hamiltonian 

^ = C \n+l,rfix'C-m,n + c\ n , n +i Uy' C m ,n + C m-l,n^l' C m,n 

+ C rn,n-l^y' C m.n + c lri,n-^ c m,n, 

U x >=a-, U y > = -a z - i\a x , F = a x . (18) 

In fact, this is the simplest scheme, exploiting non-Abelian 
hopping operators on a square lattice, which leads to a 
topological (Chern) insulating phase. Indeed, the hopping 
operators U x ^ y i should necessarily contain a term propor- 
tional to a Pauli matrix, U x \ y ' oc (y^ M , in order to induce the 
Dirac spectrum. Then, an additional term oc Xa K should be 
added to at least one of these operators in order to open the 
topologically non-trivial gap. This minimal configuration is 
realized by the operators U x / >y > in Eq. (18). 

From an experimental point of view, this scheme is signifi- 
cantly simpler to implement than the one previously exposed 
in Section III, due to the fact that it features a laser-induced 
tunneling phase in the y' direction only, and thus, it does not 
require a superlattice along x. This demands one er~ coupling 
laser that drives tunneling along x' , another laser with polar- 
ization (7 X to drive the y' tunneling, and a pair of lasers driving 
7T transitions to generate the term cx a z with the proper relative 
phase; each of this laser should also have three" sidebands" to 
drive all transitions introduced by the superlattice. This is not 
as complex as it seems at first sight, since the required side- 
bands can be generated by cascading modulators (one generat- 
ing sidebands corresponding to different Zeeman transitions, 
and another one, corresponding to the different transitions in- 
troduced by the superlattice). Nevertheless, it is fair to say 
that this represents an important technical challenge to realize 
this model experimentally. 
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(a) Pi-flux model (b) Pi-flux model with matrix hopping 




Figure 5: (a) The 7r— flux model on the square lattice. The unit cells of the n— flux lattice, with inequivalent sites A, B, are labeled by the 
coordinates (m, n) £ Z. Two inequivalent sites A, B belonging to the same unit cell are connected by a full black line. Standard nearest- 
neighbour (NN) hoppings are characterized by the tunneling factors ± J. The complex next-nearest-neighbor (NNN) hoppings, with tunneling 
factor are represented by purple dotted arrows (the hoppings with opposite factor, — iX, correspond to the reversed paths). The chirality 
introduced by the NNN terms potentially results in anomalous quantum Hall phases, (b) The same tv— flux model translated into a non-Abelian 
square lattice, with matrix hopping operators U X)V , The "undesired" diagonal hoppings D\ t i are depicted by red dotted arrows and disappear 
in the limit a — (cf. text). Note that when a — 0, this model reduces to the non-Abelian optical lattice illustrated in Fig. 1 (c). 



V. CONCLUSION 

In conclusion, we have described an experimental scheme 
to generate non-Abelian gauge potentials for cold atoms in a 
square optical lattice. The scheme generalizes the proposals 
of Refs. [29, 37] originally formulated for alkali atoms, to 
two-electron atoms (alkaline earth atoms, Ytterbium, Erbium, 
...), and requires additional superlattice potentials to generate 
finite flux per elementary cell of the lattice [43]. We have 
shown that such arrangements can be used to generate topo- 
logical phases with non-zero Chern number. We have detailed 
two examples of microscopic models exhibiting such phases, 
a first one based on an explicit mapping from the Haldane 
honeycomb lattice to a multicomponent model defined on the 
square lattice with non-Abelian tunneling operators, and a 
second one built directly on the square lattice (the so-called 
"7r-flux" model). The latter model is slightly simpler, but its 
experimental implementation would still represent a consider- 
able challenge due to the relatively large number of laser fre- 
quencies that need to be controlled precisely. As explained be- 
fore, the different transition frequencies can be generated from 
a single laser using frequency modulators, a standard tool in 
optoelectronics. As a result, we believe that the experimen- 
tal implementation is challenging but possible. Besides the 
models described here, generating non-Abelian gauge fields 



on a lattice, can simulate many different systems, including 
lattice-gauge-theory models [37] and spin-orbit coupled gases 
[7], with complete freedom over the choice of the microscopic 
Hamiltonian [33, 40, 41]. This is in contrast with Raman- 
coupled bulk systems, for example, where the form of the 
spin-orbit coupling is constrained [53]. We also stress that al- 
though the laser-coupled honeycomb lattice described in Refs. 
[23, 24] should be easier to implement experimentally, it does 
not offer the versatility of the non-Abelian optical lattices dis- 
cussed in this work. As a result, the system described in this 
paper has the potential to realize other kinds of topological in- 
sulating phases, such as, for instance, the Z2 topological insu- 
lators characterized by a topological invariant vz 2 = ±1 and 
exhibiting the quantum spin Hall effect [7, 33-36, 54-56]. 
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